Cockcroft-05-04 


Asymptotic Analysis of Ultra-relativistic Charge 

David A. Burton, Jonathan Gratus and Robin W. Tucker 

Department of Physics, Lancaster University 
and The Cockcroft Institute 

July 24, 2006 


Abstract 

This article offers a new approach for analysing the dynamic behaviour of distributions of 
charged particles in an electromagnetic field. After discussing the limitations inherent in the 
Lorentz-Dirac equation for a single point particle a simple model is proposed for a charged 
continuum interacting self-consistently with the Maxwell field in vacuo. The model is developed 
using intrinsic tensor field theory and exploits to the full the symmetry and light-cone struc¬ 
ture of Minkowski spacetime. This permits the construction of a regular stress-energy tensor 
whose vanishing divergence determines a system of non-linear partial differential equations for 
the velocity and self-fields of accelerated charge. Within this covariant framework a particular 
perturbation scheme is motivated by an exact class of solutions to this system describing the 
evolution of a charged fluid under the combined effects of both self and external electromag¬ 
netic fields. The scheme yields an asymptotic approximation in terms of inhomogeneous linear 
equations for the self-consistent Maxwell field, charge current and time-like velocity field of the 
charged fluid and is defined as an ultra-relativistic configuration. To facilitate comparisons with 
existing accounts of beam dynamics an appendix translates the tensor formulation of the pertur¬ 
bation scheme into the language involving electric and magnetic fields observed in a laboratory 
(inertial) frame. 


1 Introduction 

The intense international activity involved in probing the structure of matter on all scales, with 
particle beams and radiation, owes much to recent advances in accelerator science and technology. 
Developments in the production of high power laser radiation also offer new avenues for accelerator 
design and new diagnostic tools of relevance to medical science, engineering and the communications 
industry. A common theme in these developments is the interaction between charged particles and 
the electromagnetic field in domains where relativistic effects cannot be ignored. 

It is remarkable that many of the challenges that must be addressed in order to develop and control 
devices that accelerate charged particles have their origin in the interaction of particles with their own 
electromagnetic field. Despite the fact that the classical laws of electromagnetism were essentially 
formulated over a century and a half ago the subject of electromagnetic interactions with matter 
remains incomplete. This incompleteness has had concomitant effects on the development of quantum 
electrodynamics and renormalisation theory. At root, the difficulties reside in the recognition that 
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the quantum structure of matter at some scale is beyond observation. Furthermore, the classical 
description of the electron as a point particle leads to singularities in the Maxwell self-fields that 
inevitably create ambiguities in its interaction with the Maxwell held. The general consensus is 
that a useful domain of validity of the Lorentz-Dirac equation [1, 2, 3], describing covariantly the 
radiation reaction on a point electron, can be accommodated by performing a “reduction of order” 
that effectively replaces the equation by a perturbative second order system for the particle world 
line. One must then decide whether higher order terms in this expansion should be maintained given 
the neglect of terms associated with the regularisation scheme. This approach is behind many of the 
successful applications of approximate radiation reaction dynamics, despite the somewhat delicate 
and unsatisfactory nature of the arguments that purport to support this approach. 

It appears that the analysis of systems involving radiating matter has progressed by a very 
successful symbiosis between experimental expediency and a variety of approximation schemes with 
proven effectiveness in different domains of validity. However, as schemes for accelerating charged 
particles in these devices become more complex and ambitious in their aims it is apparent that some 
existing theoretical models are inadequate for a proper understanding of new challenges. Such models 
are sometimes enshrined in large commercial computer codes that do not survive scientific scrutiny 
or have documentation that makes contact with established scientific literature difficult. As particle 
energies increase with the use of higher intensity laser fields this issue becomes critical and more 
reliable methods for accommodating radiation reaction must be found. 

A mathematically coherent formulation of a closed system of partial differential equations describ¬ 
ing the relativistic behaviour of charged matter with electromagnetic fields is inevitably non-linear 
and, in general, exact solutions satisfying causal boundary conditions are intractable. The fact that 
the effects of (coherent and incoherent) radiation on micrometre-sized charged bunches in a host 
of newly proposed advanced devices is a significant barrier to development indicates that existing 
approaches used to model such effects are inadequate and that a new look at the whole problem is 
timely. 

Many of the radiation problems alluded to above can be circumvented by working entirely in 
the language of fields rather than both particles and fields. This makes it possible to use intrin¬ 
sic tensor held theory and exploit to the full the power of differential geometry in a relativistic 
spacetime framework. The geometric held formulation offers a number of powerful computational 
advantages over existing formulations of beam dynamics. These include the maintenance of relativis¬ 
tic covariance (local Lorentz transformation between local (accelerating) frames are unnecessary), 
the use of curvilinear coordinates adapted to the geometry of particular problems (thereby facilitat¬ 
ing imposition of boundary conditions), exploitation of available symmetries and the formulation of 
coordinate-independent approximation schemes. 

In this article a new approach is offered for analysing the behaviour of distributions of charged 
particles in a coupled electromagnetic held environment. Guided by a simple model, an approxima¬ 
tion scheme designed for charged distributions containing ultra-relativistic particles 1 is explored. The 
exact equations of motion are derived from the vanishing divergence of a relativistic stress-energy 
tensor for matter and radiation. Before presenting the approximation scheme in detail we recall the 
genesis of the radiation reaction for a point charge, mention its shortcomings and argue that these 
may be overcome by treating a large collection of charged particles as a continuum. 

1 A novel definition of an “ultra-relativistic vector field” relying entirely on the light-cone structure of spacetime is 
given in Appendix A.l. 
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1.1 Equations of motion for radiating charges 

In this paper the language of differential geometry is employed since it offers the most succinct 
mathematical framework for electrodynamics. Thus, all fields will be regarded as sections of tensor 
bundles over appropriate domains Ai of Minkowski spacetime endowed with a fixed metric g and a 
torsion-free metric compatible connection V. Sections of the tangent bundle over M. will be denoted 
YTM. while sections of the bundle of exterior p-forms will be denoted YhPM.. For any vector field 
X denote by X the associated 1-form defined by X = g(X,—). The operator d will denote the 
exterior derivative and i x the contraction operator with respect to X. To facilitate comparisons 
with existing accounts of beam dynamics Appendix A.2 translates the tensor formulation of the 
perturbation scheme into the 3-vector language involving electric and magnetic fields observed in a 
laboratory (inertial) frame. 

Due to its high symmetry, Minkowski spacetime admits a class of global charts that play a 
fundamental role in the following. A generic system of coordinates in one of these charts will be 
denoted {y /J } e M 4 where /i ranges from 0 to 3. In these coordinates the metric tensor takes the form 

g = —dy° <S> dy° + dy 1 ® dy 1 + dy 2 <g) dy 2 + dy 3 ® dy 3 . (1) 

With y° > 0 the held Jo is defined to be future-pointing. When convenient y° = ct , y 1 = x, y 2 = 
y,y 3 = z is written and the vacuum speed of light c is set to 1. The significance of this class of 
coordinates is that it offers a basis of symmetry generators {A), = J^}: 

C K ,g = o (2) 

where C x is the Lie derivative with respect to X. Such Killing vectors will be used to define energy- 
momentum densities and power-momentum fluxes associated with different held conhgurations on 
spacetime. 

If U is a domain of spacetime with boundary 

8U = + E 2 + n (3) 

for space-like hypersurfaces £i, E 2 and J a closed regular 3-form on U (i.e. dj = 0) then 

[ J = [ dj = 0. (4) 

JdU JU 

Thus 

/ J= f J- [ J- (5) 

J Si J-t .2 J n 

Such closed 3-forms imply conservation laws. 

In a vacuum the electric and magnetic helds are encoded into the 2-form T and the electric 
charge density and current are described by a source current 3-form j. The Maxwell held system on 
spacetime is 


dT = 0 

(6) 

d*g = -j 

(7) 


where * is the Hodge map associated with the spacetime metric tensor g and Q = where e 0 is 
the permittivity of the vacuum. 
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For any vector field W on spacetime and any Maxwell solution T define the “electromagnetic 
drive” 3-form 

t£? m) = X -{i w T a*T- i w *T AT}. (8) 

If W is a Killing vector field then is called a “Killing current” and 

dr$ M) = — i w T A j (9) 

eo 

follows from (6), (7) and (8). For each Killing vector field these equations establish a “local conser¬ 
vation equation” ( drw = 0) in a source free region (j = 0). 

If U is any time-like vector field (with g(U, U ) = — c 2 ) one may write uniquely 

F = eAU + *(bAU) (10) 


where g(e, U) — g(b, U) = 0 and 

^4 em) = -e A b A U - ^{g(e, e) + g(b, b)}i v (*\). (11) 

Poynting 2 related the form e A b, in a source free region, to the local field energy transmitted 
normally across unit area per second (spatial held energy current or held power) and |{#(e, e)+g(b, 6)} 
to the local electromagnetic held energy density. 

More precisely \ f E is the held energy associated with the space-like 3-chain E and ^ f s iu r u‘ ] 
is the power hux across an oriented space-like 2-chain S. 

If X is a space-like Killing vector generating space-like translations along open integral curves 
then with the split: 

4 EM) =/ixA (12) 

where iufJ>x = 0 and iuQx = 0, the Maxwell stress 2-form /j, x may be used to identify mechanical 
forces produced by a how of held momentum with density 3-form Q x ■ In any local frame {X a e 

TTM) with dual coframe [e b e TA 1 Al} the 16 functions T ( | EM ' ) = i Xb * t^ M \ where the frame 

indices a,b = 0,1,2, 3, may be used to construct the second-rank stress-energy tensor 

T (EM) = T (EM) e a 0 e b 

A standard application of these general notions is the derivation of the equation of motion of 
a charged fluid in an external electromagnetic held from the total stress-energy of the fluid and 
electromagnetic held. A thermodynamically inert (cold) fluid can be modelled with the stress-energy 
tensor 

T d) = n ^Nv®V (13) 

ce 0 

where M is a scalar number density held, m 0 some constant with the dimensions of mass, V the unit 
time-like 4-velocity held of the fluid and g(V, V) — —1. Such a stress-energy tensor gives rise to a 
set of Killing currents 

41 = — s(v,*y*(jvi6 (i4) 

Ct 0 

2 In a pre-relativistic context. 
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which are added to t ( ^ AV> to yield the total set of Killing currents for the interacting system. If one 
assumes that the electric current 3-form is j = q 0 J\f*V for some electric charge constant q 0 and that 
J\f is regular then the conservation laws 


yield the field equation of motion 


dj = 0 

(15) 


(16) 

V v V = —A V T. 
m 0 c z 

(17) 


This equation must be solved consistently with the Maxwell equations (6), (7) to determine V, Af 
and T for prescribed initial and boundary conditions. 

If the source of T contains a point charge in an external held this approach must take into account 
that part of the electromagnetic held becomes singular on the world-line of the particle so the domain 
U above must exclude this line. Furthermore, one needs to postulate a stress-energy tensor for a stable 
point particle. Without further information application of Stokes’ theorem then leaves an unknown 
contribution to the conservation law. However, experiment indicates that accelerating electrons 
(thought to be point particles) experience a reaction to their emitted radiation and that the power 
radiated is in good agreement with the covariant Larmor formula [4]. A number of methods have been 
devised in an attempt to accommodate this observed radiation reaction by adopting a regularisation 
procedure that permits the use of Stokes’ theorem in a spacetime region that includes the particle 
world-line. Inevitably this requires some assumption of how the singular Coulombic stresses are 
compensated by stresses that are not of electromagnetic origin. Early approaches invoked delicate 
limiting processes and the use of both advanced as well as retarded solutions to Maxwell’s equations 
[5, 3]. Dirac [5] offered one of the simplest covariant regularisation schemes and the resulting equation 
of motion is known as the Lorentz-Dirac equation. The use of both advanced and retarded solutions 
is however unnecessary [6, 7]. 

To highlight the assumptions in the single particle approach to radiation reaction an outline is 
given of how the Lorentz-Dirac equation can be derived using a coordinate system adapted to the 
light-cone structure of spacetime and the time-like world-line of a particle on a parametrised history 
given in the above chart as y 1 ' = £ /J (tt). For each value of u the unit tangent to the world-line is the 
4-velocity 

r = (is) 


where V ,J -{u) — ■ 

Assign to any event in spacetime not on the particle world-line, whose backward light-cone 
intersects this world-line at the event with parameter u, the set of coordinates {u, r, 6 ,0} with 
0 < r < oo,0 < 6 < n,0 < (j) < 2n using the transformation 


where for any y° > 0 : 



(19) 

L =-5jjs + t*»>+)£ = v£i 

(20) 

£>*(0,0)» = 1 

(21) 


fc=i 
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V = -g(V,L), g (V,V) = - 1, g(L,L) = 0. 

(22) 

In these coordinates the particle world-line is the curve where 3 r = 0, ( 9 , <f>) parametrise the unit 
2-sphere S 2 where {n k (9, (f))} are direction cosines and the metric tensor can be written: 

g = —e° <S> e° + e 1 <g) e 1 + e 2 ® e 2 + e 3 ® e 3 

(23) 

where 


e° = Bdu — ^ dr 

(24) 

, dr 

6 =13 

(25) 

e 2 = ^ d9 

V 

(26) 

o r 

e = — sin 9 d(j) 

(27) 

B=^l + 2r^. 

(28) 

Here, and in the following, denote d u f by / for any scalar field / and let T = C a T for any tensor 
field r. 

Suppose the total system (electromagnetic field and charge) has a stress-energy tensor giving rise 
to a regular Killing current t k such that 

dr K = 0 (29) 

in a domain U bounded by the light-like 3-chains E(n = u 0 + A u), E (u = 
3-chain E(r = ro) for arbitrary positive constants uq, Am, tq. Then 

: uq) and the time-like 

/ r K - t k + t k = 0. 

J S (u=uo +Aw) J E {u=uq ) J E (r=ro) 

(30) 

In the limit A u —> 0 this may be written 


Pjt\u 0 ,r 0 ) du + P^\u 0 ,r 0 ) du = 0 

(31) 

where 

P { K\ u o,ro)= [ f T K (u 0 ,r,9,<f>) 

Jo Js 2 

(32) 

pP(u 0 ,r 0 )= / id u T K (u 0 ,r 0 ,9,<f)) 

Js 2 

(33) 

With K e {^} and = 0 one has a similar balance of rates for each translational Killing 

vector in the basis. In terms of the “body” rate 

p( B ) _ p(B) dy „ 
dyV 

(34) 

and the “contact” rate 

p(C) ^ p(C) dyli 

dyV 

(35) 

3 For V > 0 with V time-like and g|o future-pointing the set r > 0 lies on the forward light-cone of the event ^‘ (u). 
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there is a balance of 1-forms: 

P (B) K,r 0 ) + P (C) (« 0 ,r 0 ) = 0. 

Next assume that the contact rate includes a part from the Killing currents i.e. 

p(ipI 


(36) 

(37) 


and furthermore that this part can be calculated from the retarded Lienard-Wiechert solution T to 
Maxwell’s equations for an arbitrarily moving point charge qo in no external electromagnetic held. 
In these coordinates the solution is T = dA where the 1-form 


A . go V(u) 

47re 0 r 

Integration over the 2-chain S 2 with r = r 0 yields for (33) : 


(38) 




(39) 


where A = V ,J {u)^A is the acceleration held. Clearly this 1-form is singular on the world-line where 
r 0 = 0. The hrst term however correctly accounts for the observed Larmor radiation rate of energy- 
momentum from an accelerating charge and is independent of r 0 . To proceed one must cancel the 
singular rate in P-j^ from singular terms in the remaining rates in such a way that in the limit r 0 —> 0 
the resulting system of ordinary differential equations for arbitrary uq makes both mathematical and 
physical sense. One approach is to suppose that the remaining rates are determined in terms of 
scalar helds a(u,r),/3(u,r) and the vectors V, A such that 

P (b) (m, r 0 ) + Pme C h(' u ) r o) = a(u,r 0 )V(u,r 0 ) + P(u,r 0 )A(u,r 0 ). (40) 


Inserting this in (36) and applying the projection operator II y = 1 + V A iy one hnds 

2 


^ 3 4tt e 2 + A(A} 


With this value (36) yields 

A 


d * ) go 2 1 
du l A(A) I 47reg 2r 0 


-- v \ i ‘ + lOk A(A >\ + 


2 % , 

3 4tt eg Jp) 


(41) 


(42) 


where A = V ,l {u)-^A. One way to simplify this and cancel the exposed electromagnetic singularity 
is to further assume that a(u, r 0 ) = rn ^— for some constant m 0 so that in the limit r 0 —»• 0 

one has: 


m 0 c 2 A = n V A 

3 47Te 0 

m 0 c 2 A = -x-p 5 - iv(V A A). 
o 47T6 q 


(43) 

(44) 
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This is a manifestly covariant equation of motion even though the velocity and acceleration are 
specified in a basis of Killing vector fields. Since it is tensorial it is independent of basis and in terms 
of covariant differentiation along the world-line it may be written 

m 0 c 2 A = A VvA). (45) 

O 47T6q 

This is a system of third order differential equations for the worldline £ p (u). Regarded as an initial 
value problem it requires unfamiliar initial data (£ M (0), £^(0), £ M (0)) and solutions exist corresponding 
to self acceleration which must be regarded as unphysical. 

In the presence of an external Maxwell field J- ext one must confine the motion of the particle to a 
domain U that excludes the sources of JR-xt . In this situation one expects that the equation of motion 
will acquire a contribution from the Lorentz force qoiv^ext- 

2 a 2 ~ 

m 0 c 2 A = qoivfext + x n vA (46) 

3 47T60 

so that mo is identified with the rest mass of the point particle. 

Although solutions to this system that self-accelerate can be eliminated by demanding contrived 
data at different points along the world-line there remain solutions that pre-accelerate in situa¬ 
tions where the external field is piecewise defined in spacetime [3]. Although the duration of pre¬ 
acceleration is probably classically unobservable for most electrons this general feature suggests that 
not all of the assumptions above are acceptable. 

One resolution of these difficulties [1] is to assume that the right hand side of (46) should be 
expanded as a series in q 0 with leading term for A given by t- Then the above assumptions 

are held to be accurate only to some order in q 0 and (46) is understood as: 

A —jvTex, ~ Q 2 2 :/° iv(V A VyAext) + ... (47) 

m 0 c 2 3 m 0 c 2 47re 0 

where A ex t = -^siv^ex t- The system is now manifestly a second order system of evolution equations. 
Although this offers a workable scheme it is unclear what its limitations are in different types of 
external field. Furthermore in situations where one has to contemplate the radiation from a large 
number of accelerating high-energy particles in close proximity the neglect of higher order terms in 
the expansion may be suspect. 

Given the complexities and reservations associated with (47) compared with the model leading 
to (17) the latter is adopted in this article as a description of a collection of potentially radiating 
particles in the presence of an external electromagnetic field. Thus, our fundamental set of equations 
is 


dT = 0, 

(48) 

d*F= 

(49) 

eo 

VyK= —jv?, 

(50) 

m 0 c z 

V -V = -1 

(51) 


where X ■ Y = g(X, Y ) for all vector fields X, Y e FTA4. 



To simplify the subsequent analysis the system (48-51) is recast using the 2-form F = and 

the scalar held p = where q 0 J\f is the proper charge density since M is the proper number 

density. This yields the system of held equations 


dF = 0, 

(52) 

d* F = —p*V i 

(53) 

VvV = i v F , 

(54) 

V-V= -1 

(55) 


for the triple ( V,p,F ). Equation (53) leads immediately to the integrability condition (conservation 
of electric charge) 


d*(pV) = 0. (56) 

The scalar held p will be referred to as the reduced proper charge density. Since J\f is a number 
density, J\T > 0 and so p > 0 for q 0 > 0 and q 0 < 0. 

2 Symmetric solutions to the coupled field system 

Although the held system (52-55) appears considerably less complicated than a large number of 
ordinary differential equations for a collection of accelerating charges, it is non-linear in (V, p, F) 
and so obtaining exact solutions is difficult except in certain simple circumstances. One method of 
obtaining solutions in more general situations is to employ a perturbation scheme based on properties 
of a particular class of exact solutions. 

In this section a family of exact dynamical solutions with symmetries generated by ^ and ^ will 
be described. Such highly symmetric solutions can be interpreted as accelerating “walls of charge” 
moving in vacuo under the inhuence of their self-fields (their “space-charge” helds) and an externally 
applied held. The coordinate y° will be identihed with time t in an inertial laboratory frame (the 
speed of light c = 1) and x = y\ y = y 2 , and z = y 3 . The solutions will be independent of y 1 and y 2 
i.e. 

C e F = C a F = 0 , 

dy 1 ~dy 1 

C^V = C_^V = 0 . 

Thus one may reduce the problem to a held theory on a 2-dimensional Lorentzian spacetime with 
global coordinates y°,y 3 . The held system (52-55) is solved exactly using a co-moving coordinate 
system (r, <r) adapted to the charged continuum. However, expressing the solutions in terms of 
laboratory coordinates (t,z) requires the inverse of the mapping (r, a) —> (t,z), which is generally 
difficult to obtain in closed form. To progress, a running parameter e > 0 is introduced into the 
mapping (r, <r) —> (t, z) and a perturbation scheme in e is established. This approach facilitates 
an order-by-order construction of the inverse of the mapping (r, a) —> (f, z) and thereby leads to 
1-parameter families (V £ , p £ , F £ ) of solutions in e. It will be shown that 

F‘= jr, E n F n , V = jr e"V„, ff = f >>„ 
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over some range of £ where the coefficients F n , V n and p n are 2-forms, vector fields and scalar fields 
respectively. 

Exact solutions to the system of equations (52-55) are obtained using the ansatz 


F = £(t, z) dt A dz, 

1 


yjl- p 2 {t,z) 


(d t + fj,(t,z)d z ) 


where p is the magnitude of the Newtonian velocity field of the charged continuum measured by the 
inertial (laboratory) observer d t . Equations (52) and (55) are manifestly satisfied while (53) and (54) 
lead to 


dS = p#V, _ (57) 

V y U = £#V (58) 

where # is the Hodge map associated with the volume 2-form #1 = dtAdz. Since = #(a Ah) 

for all forms a independent of dx and dip the action of iy on (57) yields 


VS = 0 


(59) 


i.e. S is constant along the integral curves of V. Let C a be a 1-parameter family of proper-time- 
parametrised integral curves of V where each value of a corresponds to an integral curve of V: 


► (t = t(r, a), z = z(t , <r)) 


v = c„e T = 8 A Si A 

OT OT 


where C' CT * is the push-forward (tangent) map associated with C a . Equation (59) is 


(C^d T )S = d T (C* a S) 
= 0 


where C* is the pull-back map associated with C a . Hence C*£ depends only on a and can be 
expressed using the function £ : M —> M where 


cm = c;s 

= £ ( t(r, a), l(r. a)) . 


(60) 


Using Vd t = Vd z = 0 the (t, z) components of (58) are 


an, , 

., . d z . . 



an, , 

.dt . . 

s^ (r ’d 

-c (<r)fr(r,<T) 
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whose particular solution satisfying the initial conditions 


t( 0, a) = 0, £(0, a) = a, 

|( 0 ,.) = 1 , |( 0 ,.) = 0 , 

at r = 0 is 

t = £(r, a) = sinh(C(<j)T), (62) 

? = o-) = ^y [cosh(C(cr)T) -1}+a (63) 

where the function ( : R —> R must be supplied as data (i.e. ( is an ingredient in the initial 
conditions). Note that the charged continuum is taken to be initially at rest in the laboratory, i.e. 

n.o= (c-a) u= a - 

More general initial conditions could have been used, but the above scenario is sufficient. 

After £ has been specified, expressions for V and F in the (t, z) coordinate system are obtained 
by inverting (62) and (63) to give (r, a) in terms of (t,z): 

t - r(t, z), 
a = a(t, z) 

and using (60) it follows that 

S(t,z) =((a(t,z)). 

An initial condition on £ is expressed using (. Equation (62) is used to show that the local 
hypersurfaces t = 0 and r = 0 are equal and so (63) yields 

<j(0, z) = z. 

Thus, (60) leads to 

£(0 ,z) = C (<r(0, *)) 

= <(*) 

on the space-like hypersurface t = 0. 

Not all choices for £ are admissible because, as noted in the derivation of (52-55), the reduced 
proper charge density p is positive: 

p = -i v (#d£) > 0 (64) 

using (57) and (55). At t = r = 0 

p(0,z) = — id t #(d t £(0, z)dt + d z £(0, z)dz) 

= d z £(0,z) (65) 

= :C'0) 

where = —dz and #dz = —dt have been used and (,'{z) = ^(z). Hence, employing (64) 

CM > 0 


11 



i.e. C is a monotonically increasing function. 

The coordinate chart (t, a) has important physical significance because it labels the streamlines 
of the charged continuum (the integral curves of V). Physically, the coordinate transformation 
(t, a) —> (t, z) is expected to be valid over the entire (t, a) plane. This will now be shown. 

The determinant A of the Jacobian matrix of the coordinate transformation (t, a) —> (t, z) is 


CV) 

C (^) 2 


+ cosh(C((j)r) ( 1 


CV) 

C (^) 2 


Since > 0 and cosh(C(cr)r) > 1 it follows that 


( 66 ) 


A > 1. 


Monotonicity of C implies that Q must vanish somewhere. Let a 0 £ M be the point where C( cr o) = 0 
and note 

lim A = 1 + Jc'Mt 2 

<T—>(T0 2 

using (66). Therefore, A is finite for all (r, a) e M 2 and the coordinate transformation (t, a) (t, z) 
is valid over the entire (t, a) plane. 

The implicit equation 

a = Z ~ + t 2 ~ X ) ( 67 ) 

for a = a(t, z) arises by combining (62) with (63). 

The initial electric field £(0, z) = ((z) = Cext + Cseif ( 2 ) where Cext models an externally applied 
field and Cseif is the self-field of the charged continuum. The field Cext is constant because it satisfies 
the vacuum Maxwell equation 

G^Cext = 0 

inside the charged continuum (consider equation (57) with p = 0). The self-field Cseif is to be defined 
in terms of the initial reduced proper charge density p(0,z): 

Cseif (■2') ~\^j P( Q ’ ^ dS ~ j z P ( 0, S ^ dS \ ( 68 ) 

and satisfies the asymptotic conditions 

Jim Cseif(2) = ^Q, 

^Jselfiz) =-\Q 

where Q = p(0, z)dz is the total charge per unit area of the continuum. 

Illustrations of the above results are shown in figures 1-4. In figures 1 and 2 the external field is 
set to zero so the charge distribution p evolves entirely under its own self-field. In figures 3 and 4 
the external field is set to a finite value. Figures 2 and 4 show the streamlines of V in spacetime for 
each case. The initial charge distribution is taken as the Gaussian 

p( 0, z) = exp(— z 2 ) (69) 
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Figure 1: The time history of the initial charge distribution (69) with zero applied external held. 


and, using (68), the associated electric self-held is 

Cseif(^) = -y^erf(z). 

The external held ("ext is chosen as: 

C “* = {0.875. 

In the hrst case (hgures 1 and 2) the initial Gaussian “wall of charge” evolves into two mirror-image 
“walls of charge” that propagate with equal Newtonian speed in opposite directions and tend towards 
the speed of light as t —> oo. In the second case (hgures 3 and 4) the initial Gaussian “wall of charge” 
is accelerated by a combination of the external held and its self-held and approximately maintains 
its shape. 

2.1 Laboratory frame description 

Although the solutions above are exact in terms of (r, a), solving (67) for a(t,z) in closed form is 
impossible if ( is arbitrary. One way of tackling this problem is to introduce a running parameter 
£ > 0, replace a by a 1-parameter family of functions a £ and solve (67) for a £ order-by-order in e. 

To effect the coordinate transformation from (r, a) to (t, z) £ is replaced with the 1-parameter 
family of functions Q £ : 

CM = J C-i + CoW- (70) 
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t 


-6 -4 -2 2 4 6 


z 

Figure 2: The lines of constant r and a generated by the initial charge distribution (69) with zero 
applied external field. Lines of constant r, except r = 0, are dark grey and lines of constant a (the 
streamlines of the flow), except a = 0, are light grey. The black line r = 0 coincides with the z-axis 
and the black line a = 0 coincides with the t- axis. 



Figure 3: The time history of the initial charge distribution (69) with finite applied external field. 
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Figure 4: The lines of constant r and a generated by the initial charge distribution (69) with finite 
applied external held. Lines of constant r, except r = 0, are dark grey and lines of constant a (the 
streamlines of the how), except a = 0, are light grey. The black line r = 0 coincides with the z-axis 
and the black curve lying in the region z > 0 is a = 0. 


where C| x t — “C-i 7^ 0 and Cseif = Co- One may interpret (70) by saying that external held effects (the 
hrst term ^C-i) dominate the effects of space-charge (the second term Co)- The coefficients of the £ 
terms in (70) have been labelled by the s-order of the term in which they appear; this convention 
will be used consistently throughout this article. 

A solution to (67) can be obtained order-by-order in £ using the series ansatz for t > 0 


= (71) 

n=0 

The electric held component £ £ {t , z) and 4-velocity held V £ can be written entirely in terms oi 
fr £ , t and z: 


£ £ {t,z ) = C £ (v £ (t,z)) 

V £ = sj\ + [C £ (a £ (t, z))] 2 t 2 d t + C £ (d £ (t, z))td z . 

Using (70), (67), (71), the choice 4 C-i > 0 and t > 0, by equating orders in e it follows that 


<?'(«, z) = z -1 + - 1 + 0(e 3 ), 


£‘(t, z) = -C-i + Co(z - t) + U + 0(e 2 ), 


(72) 

(73) 


U e — -C-i + Co(^ — t) t (d t + d z ) + 


1+2 + cs(*-o * 8 1 e+0(£2) (74) 

2^C-i C-1 


‘Similar solutions are obtained for C-i < 0 but with z + t replacing z — t. 






where Co( z ) = ^(z). Although V £ -V £ = — 1 the (t, z ) components of V £ diverge in the limits t —> oo 
and t —> 0. Divergences as t —> oo are expected because the dominant term in the electric held is 
constant, and so the charged continuum is undergoing constant acceleration to leading order in e. 
Although V £ diverges as t — * oo, the electric current J £ = p £ V £ is 

r = P £ v £ = 

= <3(* - i] (A + A) + (a , + + 0 ( £ 2 ), 

where (57) has been used. Since J £ = p £ V £ , J £ = O(e 0 ) and V £ — 0(£ -1 ) it follows that the reduced 
proper charge density p £ — 0(e). The leading order contribution to p £ is obtained by comparing (75) 
with (74) and using J 0 = piV-i where 

r=±sj„ v c = jr e n v n , f f = jr i (76) 

n =0 n=—1 n= 1 

It follows that 5 

/ = £ C^*) +0(£2) (77) 

i-lt 

Equation (75) shows that J £ is bounded to 0(e 2 ) as t —> oo. The divergence as t —> 0 of the 
coefficient of e in equation (74) stems from the square-root term in (67). Evidently, the domain on 
which the approximate solution is valid (the half-plane t > 0) is a subset of the domain of the exact 
solution (the entire (t, z) plane). This is not too surprising since the charged continuum is at rest at 
t = 0 and the leading term J 0 in the series for J £ is light-like: 


Jo ■ Jo = 0 . 


There is no observer frame in which J 0 is instantaneously at rest. 

This simple example suggests that the series 

F‘ = £ v ‘ = E £ ” v '- p’ = f >>» 

ra=—1 n =-1 n =1 

where F_i is an external held (a solution to the source-free Maxwell equations), should be inserted 
into the held system (52-55) which is then solved order-by-order in e. 

This approach can lead to the series expansions for the highly symmetric solutions just discussed, 
but can also yield £ expansions for solutions that depend on (x,y) as well as (t,z). A consequence 
of the above series expansions is that the electric 4-current has the form 

r = P £ v £ = ^Vj n . 

n =0 

It will be shown below that the above expansions partially decouple the held system (52-55) yielding 
an inhnite hierarchy of equations that are amenable to solution when supplemented with appropriate 
boundary conditions and initial data. 

5 It should be noted that the equation (f = —J e ■ V s cannot be used to calculate (77) without first calculating J 2 
since p\ = -(J 2 • V-\ + Ji ■ V 0 + Jo ■ Vi). 
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3 The perturbation scheme 

Based on the above exact symmetric solution it is assumed that there exists a class of solutions 
to (52-55) that represents configurations of charged particles in ultra-relativistic collective motion. 
It is shown in Appendix A.l that the notion of an ultra-relativistic 4-velocity vector held can be 
made precise and should be distinguished from the magnitude of a relative Newtonian speed that 
can be defined for any two time-like 4-vector fields. An ultra-relativistic velocity held is a pointwise 
limiting concept that depends on the existence of the forward light-cone structure at each event 
in spacetime. To leading order in the expansion dehned below the velocity held of the charged 
continuum is light-like. The full series will be considered as an asymptotic expansion for a solution 
to (52-55) and physically represents an ultra-relativistic conhguration controlled by the parameter 
e. Such conhgurations cannot of course be exhaustive. They are chosen to be representative of the 
class relevant to charged beams in high-energy accelerators. 

Introduce perturbation series for ( V £ , p e , F £ ) in £ of the form 


v"= 


' = E s " F » 


(78) 


where 


Theorem 1. 


V n e TTM , p n € TA 0 M , F n e TK 2 M 

(79) 

Using (52-56) the coefficients of the expansions in (78) satisfy 


dF n _! = 0 

for n G {0,1, 2,...} 

(80) 

f° 

for n = 0 


d*F n _! = l 

^ *p r Vn-r-i for n Q {1.2,...} 

(81) 

EVkJ-, 

r=0 

i = ^ ivr-iFn-r-i for n e {0,1,2,.. .} 

r=0 

(82) 

E-rV, 

r=0 

J — 1 for n = 2 

1 |0 for n G {0,1, 3,4, 5,.. .} 

(83) 

n+1 

= 0 for n G {0,1,2,...} 

(84) 


The partially decoupled equations (80-84) are amenable to an ordered analysis and can be arranged 
into the hierarchy: 


For N £ {0,1,2,...}: 

Step 3N + 1: Solve (80) and (81) with n = N for /'’ u _1 • , ^ 

Step 3AT + 2: Solve (82) and (83) with n = N for I4-i- 
Step 3N + 3: Solve (84) with n = N for p n+i . 

Proof. Equations (80-84) are obtained by inserting the series (78) into (52-55) and the integrability 
condition (56). The hierarchical structure follows by inspection. □ 
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To illustrate the partially decoupled structure in (85) consider the first 8 steps in detail: 

Step 1: Adopt an external electromagnetic held F_i i.e. a solution of the source-free equations 

dF_i = 0 and d*|C* = 0. (86) 

Step 2: Solve 

Vy jVhj — iy_ x F_i subject to V-\ ■ V_i = 0 (87) 

for l/_! ^ 0, where F_i is data obtained in the previous step. 

Although (87) is non-linear in the unknown held V_i it is straightforward to analyse; it can 
be written as a quasi-linear second-order ordinary differential equation for the integral curves 
of V-i and is the only non-linear differential equation in this scheme. 

It may seem that the pair of equations above offer 5 scalar equations for 4 unknown scalars. 
However, the contraction of the hrst equation with VLj identically vanishes, so there are only 
4 independent equations. This can be seen 6 by writing (87) with respect to an arbitrary 
basis and noting that one of the 5 components of (87) is a consequence of the other 4. 

Step 3: The leading order reduced proper charge density p\ is a solution to 

d*(piF_i)=0. (88) 

Step 4: The 2-form F 0 is a solution to Maxwell’s equations with the current piV-i as a source: 

dF 0 = 0 and d * F 0 — — * piV-%. (89) 

The conservation equation (88) ensures that (89) has a solution. 

Step 5: The zero order velocity held V 0 is obtained from the linear equations 

Vv_jVo + Vy 0 V_i = iy.jFo + iv 0 F- 1 subject to V-\ ■ Vo — 0. (90) 

Step 6: The second order coefficient p 2 in the £ expansion for the reduced proper charge density p £ 
satisfies 

d*(p 2 TAi) + d*( Pj f p ) =0. (91) 

Step 7: Equation (91) ensures that 

dFi = 0 and d * /-j — — * p 2 E_i — *piV 0 (92) 

may be solved for the coefficient F\ (the hrst order term in the £ expansion of F £ ) given 
appropriate initial and boundary conditions. 

Step 8: The linear partial differential equation for Vi is 

V v _ 1 V± + Vvgho + Vy x Vki = iy^Fx + i Vo F 0 + i Vl F_ 1 subject to 2E_i • lA + V 0 ■ V 0 = -1. 

(93) 

6 The split of (87) with respect to the laboratory frame is given in Appendix A.2. The equation obtained by 
contracting (87) with the laboratory observer field d t is redundant because it follows as a consequence of (119) and 
( 120 ). 
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V £ must be calculated to at least first order in e to obtain a time-like vector field. To see this 
let X £ — £~^V—i + V 0 and note X £ ■ X £ = V 0 ■ V 0 using the metric-product conditions in (87) and 
(90). Using lemma 5 in Appendix A.l Vo • Vo > 0 and so X £ is light-like or space-like. Now introduce 
Y £ =X £ + eV x = + U 0 + eVi and note Y £ -Y £ = V 0 -V 0 + 2Vi ■ V^ + 0(e) = -1 + 0(e) using 

(93). Therefore, there exists some value of e > 0 for which Y £ is time-like. 

It has already been stated that the pair of equations (87) in step 2 may be considered as 4 
non-linear equations in 4 unknowns. Furthermore: 

Lemma 2. The equations in each step 3N + 2, for N > 1, are equivalent to a set of 4 independent 
inhomogeneous linear equations for the 4 unknown components o/Vn-i- 

Proof. Step 377 + 2 involves solving (82) and (83) for Vn- i- Assume that the equations in steps 1 to 
31V + 1 have been solved and define L n as 

L n = (v Wl K -1 - iv^Fn-r-l) ■ (94) 

r =0 V 

Clearly (82) is equivalent to L n = 0. The vector fields U_i,..., Vn-i have been obtained in steps 1 
to 31V + 1 and so L n = 0 for n = 0,..., AT — 1. It will now be shown that given any Uv-i satisfying 
(83) with n = N the scalar field Iv-iLn identically vanishes. 

From (83) for n = 0,..., N 

»=(xa-i-v^A -£(76^,) • w, 

\r=0 ) r =0 r =0 

n 0 

= • Vn-r-l + E ■ Vr.v. 

r=0 q=n 

where the second sum has been re-labelled by q = n — r. The two sums are identical so 


• U r _i = 0. (95) 

r =0 


Using (94) and L n = 0 for all 0 < n < N 

N 

fy-l-bjv =y^fy„_i-bjy_ n 

n =0 

N N-n 

= ^fy n -! ^2 (Wjv-n-r-iK-l ~ iv r -i F N-n-r-2j 
n =0 r=0 V 

N N-n 

= (^-1 • - iVn-M-iFN-n-r-x) 

n =0 r =0 
N N 

= ( V »- 1 • ^v N -,-iV a -n-$, - iv^ivs-n-iFN-s- 1) (96) 

n =0 s=n 


where, in the final step, the innermost sum has been re-labelled by s = n + r. The double summation 
in (96) is over the set 

{( n , A') | n. .s e. {(). 1. N}, n < s} 
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and so 


N N Ns 

£X>- = £X>‘ 

for any summand S ns . Therefore, using (96) 

iv-iL N = ^ | ^14-i • Vv^K-n-i — d'„ n - l I 

s=0 \n=0 n=0 / 

= 0 (97) 

where the first term in the summand vanishes because of (95) and the second term vanishes because 
= — ivs-n-ii-Vn-i and the sum over s is from 0 to n. 

Note that the lemma applies to step 5 and beyond and so, for the purposes of the lemma, the fields 
F_i, F 0 obtained in steps 1 to 4 are prescribed data. Let {Xi, X 2 , X 3 , X 4 }, with X t e YTM. 

for i = 1, 2, 3,4, be any vector frame such that Xi = U_i and define {e 1 , e 2 , e 3 , e 4 }, where e l e , 

to be the coframe dual to {X\, X 2 , X 3 , X 4 }: 

Axd - s} 

where 5] is the Kronecker delta. Equation (82) can be written as 

C{V N - 1 ) 6v-i (98) 

where C is the R-linear differential operator on TA 1 At 

C(V N - 1 ) - V v - ,14-* + + iv^F-i 

where £n-% € VXM, depends on {U_i, V 0 ..., V N - 2 , F 0 ,..., F n _i} but is independent of Ev-i. 
Hence, using (97) 

*Xi(£(Vw-i) — £n- 1 ) = i-v , (-C(Kv-i) — 6v-i) — 0 (99) 

Likewise, equation (83) can be expressed as 

ix^N-i = iv-iVn-i = Cn- 1 ( 100 ) 

where the scalar held C N _ 1 e rA°A4 depends only on {V 0 ,..., V N _ 2 }. 

Using the coframe {e 1 , e 2 , e 3 , e 4 } 

4 

Vn-i = ^ 

for component fields Ajv-i,*- From (100) Xn-i,i = Cn- 1 and (98) becomes 



^(CW-re 1 ). 


( 101 ) 
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From (99) ix x annihilates (101) identically and 


ixj f A/y-i^-e*) J — lx, (&v- i — for j — 2,3,4. 

Thus, using the inhomogeneous linear algebraic constraint (83), 3 inhomogeneous linear differential 
equations have been obtained for the 3 unknown fields {Ajv_i )2 , Xn-i, 3 , A.v- 1 , 4 } and Ajv- 1,1 is obtained 
using (83), determining fov-i- □ 

In conclusion, the only non-linear equations in the hierarchy developed in theorem 1 are those 
for Hi. All of the remaining equations are linear and therefore, in principle, simpler to solve. 

4 Example : Charged beam propagating in free space 

As a simple application of the hierarchy (80-84) consider a high-energy charged beam propagating 
in free space. Unlike the exact “wall of charge” solutions derived earlier, the charge distributions 
considered here have finite extent transverse to and along the direction of propagation. The first 
eight steps of the hierarchy (80-84) are: 

Step 1: The leading order Maxwell field F_ 1 (the external field) is chosen to vanish so (86) is trivially 
satisfied. 

Step 2: Using (87) the leading order 4-velocity field U_i satisfies 

V = 0, U_1 • U_! = 0 (102) 

since the external field F_i = 0. Therefore, consider solutions to (80-84) adapted to a null 
geodesic coordinate system (u, v, x, y) where u = z — t, v = z + t and 

g = —dt <E> dt + dx ® dx + dy ® dy + dz <g) dz 

= ^ (du <g) dv + dv (8) du) + dx <S> dx + dy <g) dy (103) 

with the spacetime volume 4-form 


*1 = dt Adz A #_l1 


where 


#_l1 = dx A dy. 


The pair (d t , d z ) and their duals (dt, dz) are related to (d v , d u ) and (dv, du) as follows: 


s„ = l(0,+a,), d u = \(d z -d t ), 

dv = dz + dt, du = dz — dt. 

Metric identities used in the subsequent steps are 

*(du A a) — (— l) p du A #_lo:, *(dv A a) = — (— l) p dv A #±a (104) 
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where #±(a AX) — ix#±ot and a is any p-form such that ig u a = ig v a = 0. The vector field 
V-i is chosen as 

V-i = 7-1 (dt + d z ) 

= ^l-id v (105) 

where 7 _i ^ 0 is assumed and, using Vd v — Vd u = 0, the only non-zero component ol (102) 
is 

d„7_i = 0. (106) 

Hence, 7 _i is independent of v. 

Step 3: Using (105), equation (88) is 

9 v (pij-i) = 0 

and so, using (106) 

d v pi = 0 . 

It follows that p\ is independent of v. 

Step 4: Using (104), particular solutions to (89) may be written 

F 0 = o A du (107) 

where the scalar potential $ 0 is chosen to satisfy 

d v $ 0 = 0, o = Pi7-i#±l (108) 

and d± is the exterior derivative in the (x, y ) plane defined by 
dj_f — d x f dx + d y f dy 

for any 0-form / <G TA°A4. 

Step 5: Using (90), the equations for V 0 are 

Vy_ 1 Vo + Vy 0 U_ 1 = iv_ 1 Fo, %! ■ Vo = 0 
and are solved by the particular solution 

U 0 = 0 

since iv^Fo = 0 using (107), (108) and (105). 

Step 6: Since V 0 = 0 it follows that (91) has the form 

d * (p 2 U-i) = 0 

and, using (105) and (106), is solved by any p 2 satisfying 
d v p 2 = 0 

i.e. p 2 is independent of v. 
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Step 7: Since V 0 = 0 the Maxwell equations (92) are 

dFi = 0, d* Fi — —p 2 * V-i 

and, employing the same method used in step 4, are solved by the particular solution 

F l = d$i A du (109) 

where is chosen to satisfy 

d v $i = 0, di <I>, ' P27 -i#±!- (HO) 

Step 8: Using (93), the equations for V\ are 

V y_i hi + Vvilti = i v ,/q. VUi • Vi =■ (111) 

and since iy^Fi — 0 (which follows using (109), (110) and (105)) a particular solution to 
(111) is 

Vi = ~ 7 ^—d u , a u7 _x = 0 

27-i 


using (106), (105) and V<9„ = V d u = 0. Hence, the function 7 _i is independent of u and 
therefore, using (106), 7 _i depends only on x and y. 

In summary these solutions, describing an ultra-relativistic charged distribution propagating along 
the z-axis with its electromagnetic self-fields, are 

V £ =-2 j^dv-e-Lm + Ole 2 ) 

£ Z 7 _i 

= (a - 1 filyj St + (a - 1 fayr)+° <£2) 

F £ = o A du + ed$i A du + 0(e 2 ) 

= — (d$ 0 + ed$i) A dt + (d$ 0 + eddq) Adz + 0(e 2 ) 

where 4> 0 and $i satisfy 

d±#±d_L$ 0 = 7 _iPi#_l 1, d_L#_Ld_L$i = 7 -ip 2 #±l (112) 


and 


$0 = x, y ), 4>i = (u, x , y), 

pi = Pi(u,x,y), p2 = p2(u,x,y), 

7-i — y) 

with u = z — t. The scalar fields p 0 ,Pi,7-i are determined by their values on the space-like hyper¬ 
surface t = 0 given as data. The potentials d> 0 and dq are then solved in terms of 7 _ipo and 7-iPi 
using the 2-dimensional Poisson equations (112) in the (x,y) plane. 

The 3-velocity of the beam is along the direction z in the laboratory frame with Newtonian speed 
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For example, consider a Gaussian bunch with transverse radius Ro travelling at constant Newtonian 
speed 1 — to order e 2 \ 

7-i (x,y) = bo, pi(z,x,y) = a 0 exp 

where ao, Ro and bo are constants and S : M —> M is a smooth bump function vanishing outside the 
interval (— Zi,Zi) and S(z) = 1 for z E (—z 2 ,z 2 ) and zi > z 2 > 0. Then the laboratory reduced 
charge density 7 _ipi for some range of t is 



/ _|_ y 2 \ 

7-iPi =l-i(x,y)pi( z ~t,x,y) = a 0 b 0 exp (- JE(z-t). 

Working in the cylindrical polar coordinates (t, R, (j), z) where x = R cos (j) and y = R sin 0, a cylin- 
drically symmetric solution to (112) well-behaved at R = 0 is 

and the corresponding electromagnetic 2-form F 0 is 

r ? r / r ? \ i 

F 0 = dobo ^ |^1 - exp JJ — t) dR A (-dt + dz). 

The laboratory electric held is radial, the magnetic held is azimuthal and their magnitudes are equal 
and vanish outside of the support of S. 

5 Conclusion 

When analysing the dynamics of charged particle beams it is often fruitful to adopt a description 
based on classical helds rather than classical point particles. Pathologies (such as pre-acceleration) 
associated with radiating point particles are avoided by relying on held-theoretical notions. 

A novel analysis of charged beam dynamics has been presented and a model of a freely propagating 
charged bunch discussed. The approach relies on an asymptotic series representation of solutions 
to self-consistent spacetime covariant held equations for (V, p, F) describing a charged continuum. 
The asymptotic series for V is based on the light-like vector held V-\ leading to an ultra-relativistic 
approximation. The hierarchy of equations obtained are more amenable to analysis than the original 
non-linear held system and particular solutions have been presented. 

There are numerous avenues for the development of this work involving ultra-relativistic charged 
beams in the vicinity of beam pipes, RF cavities, spoilers, etc. leading to dynamical effects that are 
often described in terms of “wake-helds” [8, 9]. This work will lead to a clearer understanding of 
radiation-reaction exhibited by continuum models of charged particle beams. 
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A Appendix 

A.l Definition of ultra-relativistic and light-like limited vectors 

The spacetime description of a classical point particle with mass mo > 0 invokes a parametrised 
curve (world-line) with a future-pointing time-like tangent vector. For affinely parametrised curves 
the tangent vector V (4-velocity) is normalised: V • V — — 1. A spacetime reference frame may be 
associated with a time-like vector field U with U ■ U = —1 and an “observer” in such a frame is 
modelled by an integral curve of U. The 4-momentum of the particle is defined to be p = m 0 V. 
At an event where an observer curve intersects the world-line of the particle one has the orthogonal 
decomposition: 


p = £U + P 


where U ■ P = 0. Relative to U, the scalar £ is the energy of the particle (in units with c = 1) and 
P is its 3-momentum. The Newtonian velocity v of the particle relative to U is 


v = 


P 

£ 


(113) 


Since U ■ U = V ■ V = — 1 implies £ 2 = P • P + it is traditional to say that the particle is 
“relativistic” relative to U if £ 2 ~ P • P or equivalently v • v ~ 1. Clearly, such a notion depends 
on both V and U and the nearness of hJ- to unity. Furthermore, one may always find a frame in 
which P = 0 and £ = m 0 . However, given a particular frame U one may contemplate a 4-velocity V 
that is relativistic relative to it. To generalise this notion for any frame one may consider a family 
of 4-velocities V £ and invoke properties of V s as £ tends to some limit. 
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Let £ e (0, £ max ) = {e G M|0 < s < £ ma x} be the running parameter introduced in the main body 
of the article. A 1 -parameter time-like vector field (which is not necessarily normalised) W £ is 
ultra-relativistic if for any time-like vector field Z e TTXl 


lim 

6->0 


( W £ -z | \ 

V IM ) 


= oo 


(114) 


where ||AT|| = \J—X ■ X for any time-like (or light-like) vector held X e TT XI. Note that in this 
definition Z is not a 1-parameter vector held. 

Lemma 3. Our definition of an ultra-relativistic vector field (Ilf) on spacetime is independent of 
the choice of the time-like vector field Z. 

Proof. Clearly (114) is independent of any rescaling of Z. Therefore let Z and Z be two ^-normalised 
future pointing time-like vector helds (i.e. Z ■ Z = Z ■ Z = — 1) and dehne the vector held X\ where 


Hence 


(Z+(Z-Z)Z) 
((Z ■ Z) 2 - l) 1 / 2 ' 


Z = -(Z ■ Z)Z + ((Z ■ Z) 2 - l) 1/2 Xi 


(115) 


and X\ ■ Z — 0, X% ■ X\ = 1. Choose space-like unit vector helds X 2 , A 3 so that (Z, X l} X 2 , X 3 ) is s 
r/-orthonormal frame. It follows that 


which implies 


W £ ■ W £ = -(W £ ■ Z ) 2 + {W £ ■ X r ) 2 + (W £ ■ X 2 ) 2 + (W £ ■ X 3 ) 2 


(W £ • xfi ) 2 _ (W £ • Z) 2 1 (W £ ■ X 2 ) 2 {W £ ■ X 3 ) 2 

Ill'll 2 “ ll^ll 2 IM 2 Ill'll 2 


im 


Ill'll ■ 


The metric contraction of (115) with is 


n^'ii ( ] n^'ii (( ] ] n^'ii 


(116) 


^Z-W £ Z-VT £ t/f , ^ 2 


which, using the triangle inequality for scalars, implies 


IZ < 11X11 1 ,(Z Z) 2 

1 \\W‘\\ ~ \\W-\\ +(( > > \\W'\\ 
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Rearranging the above equation gives 


\Z-W £ \ 




jZ-W £ l 

“p^ 


- (( ZZ ) 2 - l) 1 / 2 


\X x -IF £ | 

\\W'\\ 


and so, using (116) 

l|W £ ll - 1 1 \\W'\\ (( 


' z > 2 - w=W 0* • ■ 21 - ({z ■ z? - 1)1/2 )' 


Since \Z-Z\- ((Z ■ Zf - l ) 1 / 2 > 0, and Z and Z are independent of e, it follows that 


. IZ-W £ I 

“i w^r 


• lim 


\Z-W £ \ 

■« ill'll 


The connection between our definition of an ultra-relativistic vector field and v v ~ 1 is exhibited 
by introducing the 1-parameter family of Newtonian velocities v £ of W £ with respect to the observer 
field U : 


where U -U — —1. It follows that 


W 

W^-U 


W £ 

v - [ w s -U + 


W £ 

W e -U~ 


W £ W £ 2(W £ -U) TT TT W £ -W £ „ 

- 1 —i- - + u ■ U =-h 1 

(w* ■ uy w £ -u ( it £ • uy ' 


Thus lim £ ^o v £ • v £ = 1 <*=* lim £ ^ 0 W £ ■ W £ /(W £ ■ Uf = -lim £ ^ 0 (\\W £ \\ /\W £ ■ U\f = 0 and 
hence lim £ ^ 0 v £ • v £ = 1 <*=>• lim £ ^ 0 | W £ ■ C/|/||1F £ || *->■ oo (i.e. W £ is ultra-relativistic). 

Observe that in definition (114) lim £ ^ 0 W s may or may not exist. Given a 1 -parameter vector 
field, W £ , W £ is said to be light-like limited if W° = lim £ ^ 0 W £ exists and is light-like (i.e. W° ^ 0 
and 11 W° 11 = 0). The relationship between ultra-relativistic vector fields and light-like limited vector 
fields is given by the following lemma: 


Lemma 4. If W £ is light-like limited then W £ is ultra-relativistic. If W £ is ultra-relativistic and the 
limit lim £ ^ 0 W £ exists and is nowhere 0 then W £ is light-like limited. 

Proof. If W £ is light-like limited then lim £ ^ 0 W e • Z\ = |lim £ ^ 0 (W £ ) • Z\ — |W° ■ Z | > 0 since 
W° is light-like and Z is time-like. Furthermore, lim £ ^ 0 (>V £ • W £ ) = — ||W °|| 2 = 0 and so W £ is 
ultra-relativistic. 

If W £ is ultra-relativistic and the limit W° = lim £ ^ 0 W £ exists and is nowhere 0 then (since W° • Z 
is hnite) ||W°|| must vanish to ensure (114) is obeyed. □ 

The 1-parameter family of velocity fields V £ introduced in the main body of this article (see section 
3) is ultra-relativistic. Moreover, the vector fields V_i, Vq and Vi then have certain properties: 


Lemma 5. 


• The 1-parameter family of vector fields V £ is ultra-relativistic if and only ifV -1 does not vanish 
anywhere. 
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If V £ is ultra-relativistic then Vq is nowhere time-like and V\ does not vanish anywhere. 


Proof. 

• If V-i 7 ^ 0 then for any time-like vector Z, V-i ■ Z ^ 0 so 

lim(|^ • Z\/ ||F1) = linr^- 1 !^! ■ Z\) = oo 

since ||T £ || = 1 and so V £ is ultra-relativistic. If there is some x <G M. such that V_i(x) = 0 
then 


lim(|I/ £ (x) • Z(x )|/ ||T £ (x)||) = lim(£ 0 |y 0 (x) • Z{x)\) ± oo 
so V £ is not ultra-relativistic. 

• Assume V 0 is time-like. Since V £ is ultra-relativistic it follows VLi ^ 0 and, furthermore, since 
Vq is time-like VIi • V 0 ^ 0 ; however, this conclusion contradicts (90) which states V 0 and F_i 
are orthogonal. Hence, Vo is not time-like (i.e. Vo is space-like, light-like or zero). 

Now since VLi is light-like and Vo is not time-like then in order to satisfy 2V_\ -V\ +Vq-V 0 = —1 
from (93) V% does not vanish anywhere. 


□ 


A.2 The perturbation hierarchy in Gibbs’ 3-vector notation 

To facilitate a comparison with existing fluid descriptions of beam dynamics [10, 11] that employ 
Gibbs’ 3-vector notation the equations (52-56) in section 3 are here transcribed into that language. 
By splitting with respect to an inertial observer field \d t where 

g = —c 2 dt ® dt + dx ® dx + dy <g) dy + dz ® dz 


one obtains 


and 


m 0 c 

V • e =-7p, 

Qo 

__ , 1 1 de 

V xb= -PP+- —, 

Qo 

„ db 

V x e+- = °, 

V • b = 0, 


2 dt ’ 


p 


dp 

1 t , 

dt \m 0 

-f + - 


V)p = qo ( 1 e+-p> 


mo-Qt^p) + V ' (PP) • 0 


a , P I-? Qo 

idtTt +-Vy = —^p ■ 

m 0 m,QC z 


(117) 


(118) 
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where a ■ b is the standard Euclidean dot product of two 3-vector fields a and b. Note that (118) can 
be obtained by taking the scalar product of p with the differential equation for p and the algebraic 
constraint on 7 and p in (117). Therefore, like the final equation in (117) (charge conservation), 
equation (118) is redundant. However, the two differential equations are treated very differently in 
the subsequent asymptotic analysis. Equation (118) is automatically maintained order-by-order by 
the algebraic constraint on 7 and p but charge conservation must be imposed order-by-order. Hence, 
(118) is not included in the following. 

The scalar field p is the reduced proper charge density introduced in the main body of the article 
and the MKS fields (p, e, b ) and 7 are related to (V, F) by 


7 P 

V — -d t +- 

c m 0 c 

F = dt A E + * (dt A m) 

m 0 c m 0 V / 

where (P,E,B) are the vector fields 

P = (p-i)d x + (p-j)d y + (p-k) d z , 

E — (e ■ i) 0 X -i- (e ■ j) d y + (e ■ k) d z , 

B = (6 • i) d x + (b ■ j) d y + (b ■ k ) d z 

and (i,j,k) is an orthonormal inertial 3-vector frame 
The expansions for the 1-parameter families e £ , b 

e £ — e 1 + Cq + £C\ + ... , b £ - 
7 £ — £ _1 7-i + 7o + £7i + • • • , P £ 
p £ = ep x + e 2 p 2 + ... 

Substituting the above expansions into (117) yields the following hierarchy corresponding to equations 
(86-93): 

Step 1: The source-free Maxwell equations are satisfied by the prescribed external electric held e_i 
and external magnetic held b p. 


: , 7 £ , p £ and p £ corresponding to (78) are 
= £ 1 b_i + 60 + £b\ + . . . , 

= £ _ 1 p_i + po + £Pi + ■ ■ ■ 




(119) 

( 120 ) 
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Step 3: The leading order equation describing charge conservation is 

^( 7 -iPi) + V ' = 0 

and is solved for p \. 

Step 4: Solve the Maxwell equations for e 0 and 6 0 with 7 _i/?i and ^-piP-i as sources: 


m 0 c 

V • e 0 =-7-iPi, 

% 

T7 n 1 , Ideo 

V x 6 0 = — PiP-i + 

go c 2 at 

V • b 0 = 0. 


Step 5: Solve 


dp 0 dp_ i 

7 - i ^r +T "^r 


= g 0 (7-i e o + Toe-0 H- (p ~i x 6 0 + p 0 x 6 _i) 

mo 


- 7_i 70 + P \ f ° = 0 

mgC 2 

for 70 and p 0 Use of the above algebraic equation to eliminate 70 leads to an inhomogeneous 
linear partial differential equation for p 0 . 

Step 6 : The equation 

m 0 J^(7oPi) + m 0 J^(7-ip 2 ) + V • (pipo) + V • (P2P-1) = 0 
enforcing charge conservation is solved for p 2 . 

Step 7: The electric and magnetic fields e\ and 61 are solutions to the following Maxwell equations: 

m 0 c 2 


V • ei = - 


- ( 7 - 1 P 2 + ToPl) i 


Qo 

„ , 1 , x 1 de x 

Vx bl = -(p 2 p_ 1+ p lPo ) + - 2 -^, 

Vxe 1 + m.o, 

V • 61 = 0. 


Step 8: Solve 


7 -w 1 + 7 »l 2+7 '% 1 + f-' v )p‘ + f-' v )p» + f-' v )^ 

at at at \m 0 ) \m 0 J \m 0 / 


- 2 q_i 7 i - 7070 + 


= g 0 (7-iei + 7 oe 0 + 71^-1) + — (P- 1 x 61 + p 0 x 6 0 + pi x 6 _i), 
m 0 

2 P-i -Pi+Po-Po _ 
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for 71 and p\. Eliminating 71 using the above algebraic equation leads to an inhomogeneous 
linear partial differential equation for p\. 


31 



